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Abstract 

We investigate hadron production as well as transverse hadron spectra in nucleus-nucleus col- 
lisions from 2 ^-GeV to 21.3 A-TeV within two independent transport approaches (UrQMD and 
HSD) that are based on quark, diquark, string and hadronic degrees of freedom. The compar- 
ison to experimental data demonstrates that both approaches agree quite well with each other 
and with the experimental data on hadron production. The enhancement of pion production in 
central Au+Au (Pb+Pb) collisions relative to scaled pp collisions (the 'kink') is well described by 
both approaches without involving any phase transition. However, the maximum in the /ir^ 
ratio at 20 to 30 A-GeV (the 'horn') is missed by ~ 40%. A comparison to the transverse mass 
spectra from pp and C+C (or Si+Si) reactions shows the reliability of the transport models for 
light systems. For central Au-|-Au (Pb+Pb) collisions at bombarding energies above ~ 5 A-GeV, 
however, the measured m^-spectra have a larger inverse slope parameter than expected from 
the calculations. The approximately constant slope of spectra at SPS (the 'step') is not re- 
produced either. Thus the pressure generated by hadronic interactions in the transport models 
above ~ 5 A-GeV is lower than observed in the experimental data. This finding suggests that the 
additional pressure - as expected from lattice QCD calculations at finite quark chemical potential 
and temperature - might be generated by strong interactions in the early pre-hadronic/partonic 
phase of central Au+Au (Pb+Pb) collisions. 

1 Introduction 

The phase transition from partonic degrees of freedom (quarks and gluons) to interacting hadrons is 
a central topic of modern high-energy physics. In order to understand the dynamics and relevant 
scales of this transition laboratory experiments under controlled conditions are presently performed 
with ultra-relativistic nucleus-nucleus collisions. Hadronic spectra and relative hadron abundancies 
from these experiments reflect important aspects of the dynamics in the hot and dense zone formed in 
the early phase of the reaction. Furthermore, as has been proposed early by Rafelski and Miiller 
the strangeness degree of freedom might play an important role in distinguishing hadronic and partonic 
dynamics. 

In fact, estimates based on the Bjorken formula ^ for the energy density achieved in central Au+Au 
collisions suggest that the critical energy density for the formation of a quark-gluon plasma (QGP) is 
by far exceeded during a few fm/c in the initial phase of the collision at Relativistic Heavy Ion Collider 
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(RHIC) energies jS], but sufficient energy densities (~ 0.7-1 GeV/fm^ might already be achieved at 
Alternating Gradient Synchrotron (AGS) energies of ~ 10 A-GeV jniEl- More recently, lattice QCD 
calculations at finite temperature and quark chemical potential /ig j7] show a rapid increase of the 
thermodynamic pressure P with temperature above the critical temperature Tc for a phase transition 
to the QGP. The crucial question is, however, at what bombarding energies the conditions for the phase 
transition (or cross-over) might be fulfilled. 

Presently, transverse mass (or momentum) spectra of hadrons are in the center of interest. It is 
experimentally observed that the transverse mass spectra of kaons at AGS and SPS energies show a 
substantial flattening or hardening in central Au+Au collisions relative to pp interactions (cf. Refs. jHl 
in]). In order to quantify this effect, the spectra are often parametrised as: 



where uit = y 'w? + p^ is the transverse mass and T is the inverse slope parameter. This hardening of 
the spectra is commonly attributed to strong collective fiow, which is absent in the pp or pA data. 

The authors of Refs. piU ITT] have proposed to interpret the approximately constant K"^ slopes 
above ~ 30 A-GeV - the 'step' - as an indication for a phase transition following an early suggestion 
by Van Hove ^2]- This interpretation is also based on a rather sharp maximum in the j-n^ ratio at 
~ 20 to 30 A-GeV in central Au-|-Au (Pb+Pb) collisions (the 'horn' fOlE])- However, it is presently 
not clear, if the statistical model assumptions invoked in Refs. [TIH ITT] hold to be reliable. 

We will demonstrate in this contribution that the pressure needed to generate a large collective fiow 
- to explain the hard slopes of the spectra as well as the 'horn' in the j 7r+ ratio - is not produced 
in the present models by the interactions of hadrons in the expansion phase of the hadronic fireball. 
In our studies we use two independent transport models that employ hadronic and string degrees of 
freedom, i.e., UrQMD (v. 1.3) p!3t IT^ and HSD I16j . They take into account the formation and 
multiple rescattering of hadrons and thus dynamically describe the generation of pressure in the hadronic 
expansion phase. This involves also interactions of 'leading' pre-hadrons that contain a valence quark 
(antiquark) from a primary 'hard' collision (cf. Refs. [T3| I17j). 

The UrQMD transport approach 1^ includes all baryonic resonances up to masses of 2 GeV as 
well as mesonic resonances up to 1.9 GeV as tabulated by the Particle Data Group [TH]. For hadronic 
continuum excitations a string model is used with hadron formation times in the order of 1-2 fm/c 
depending on the momentum and energy of the created hadron. In the HSD approach nucleons, A's, 
N*(1440), N*(1535), A, S and S* hyperons, S's, S*'s and fi's as well as their antiparticles are included 
on the baryonic side whereas the 0~ and 1~ octet states are included in the mesonic sector. High energy 
inelastic hadron-hadron collisions in HSD are described by the FRITIOF string model whereas low 
energy hadron-hadron collisions are modeled based on experimental cross sections. Both transport 
approaches reproduce the nucleon-nucleon, meson-nucleon and meson-meson cross section data in a 
wide kinematic range. We point out, that no explicit parton-parton scattering processes (beyond the 
interactions of 'leading' quarks/diquarks) are included in the studies below contrary to the multi-phase 
transport model (AMPT) |20j, which is currently employed from upper SPS to RHIC energies. 

2 Hadron excitation functions and ratios 

2.1 "p-p versus central AA reactions — the 'kink' 

In order to explore the main physics from central AA reactions it is instructive to have a look at the 
various particle multiplicities relative to scaled pp collisions as a function of bombarding energy. For 
this aim we show in Fig. ^the total multiplicities of 7r^,i^^ and K~ (i.e., the 47r yields) from central 
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Figure 1: Total multiplicities 
of 7i'~^,K~^ and K~ (i.e., 47r 
yields) from central Au+Au 
(at AGS) or Pb+Pb (at SPS) 
collisions in comparison to the 
total multiplicities from pp 
collisions (scaled by a factor 
350/2) versus kinetic energy 
E'lab- The solid lines with full 
triangles and squares show the 
UrQMD (l.h.s.) and HSD re- 
sults (r.h.s.) for AA colli- 
sions, respectively. The dotted 
lines with open triangles and 
squares correspond to the pp 
multiplicities calculated within 
UrQMD (l.h.s.) and HSD 
(r.h.s.). The figure is taken 
from Ref. |17j . 



Au+Au (at AGS) or Pb+Pb (at SPS) collisions (from UrQMD and HSD) in comparison to the scaled 
total multiplicities from pp collisions versus the kinetic energy per particle -Eiab- 

The general trend from both transport approaches is quite similar: we observe a slight absorption 
of pions at lower bombarding energy and a relative enhancement of pion production by rescattering in 
heavy-ion collisions above ~10 A-GeV. Kaons and antikaons from AA collisions are always enhanced 
in central reactions relative to scaled pp multiplicities, which is a consequence of strong final state 
interactions. Thus, the 'kink' in the pion ratio as well as the enhancement might result from 
conventional hadronic final state interactions. 



2.2 Particle yields in central collisions of heavy nuclei 

Fig. 121 shows the excitation function of tt^ ,tt^ , , K~ and A + S° yields (midrapidity (l.h.s.) and 
rapidity integrated (r.h.s)) from central Au+Au (Pb+Pb) collisions in comparison to the experimental 
data ^. As can be seen from Fig. |21 the differences between the independent transport models are 

^Note that all data from the NA49 Collaboration at 30 A-GeV have to be considered as 'preliminary' 
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Figure 2: The excitation func- 
tion of TT^ , IT" , , K~ and 
A + E° yields from 5% cen- 
tral (AGS energies, SPS at 160 
A-GeV and at RHIC energies), 
7% central (20, 30, 40 and 80 
A-GeV), 10% central for A+S° 
at 160 A-GeV Au+Au (AGS 
and RHIC) or Pb+Pb (SPS) 
collisions in comparison to the 
experimental data from Refs. 

|nil2ai2Sj (AGS), [2111231211 

(SPS) and [211 l2Hll2i (RHIC) 
for midrapidity (left column) 
and rapidity integrated yields 
(right column). The solid lines 
with open squares show the re- 
sults from HSD whereas the 
dashed lines with open trian- 
gles indicate the UrQMD cal- 
culations. The lower theoret- 
ical errorbars at RHIC ener- 
gies correspond to the yields 
for 10% central events. The fi- 
gure is taken from Ref. 
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less than 20%. The maximum deviations between the models and the experimental data are less than 
~ 30%. In addition, a systematic analysis of the results from both models and experimental data for 
central nucleus-nucleus collisions from 2 to 160 A-GeV in Ref. |T7| has shown that also the 'longitudinal' 
rapidity distributions of protons, pions, kaons, antikaons and hyperons are quite similar in both models 
and in reasonable agreement with available data. The exception are the pion rapidity spectra at the 
highest AGS energy and lower SPS energies, which are overestimated by both models For a more 
detailed comparison of HSD and UrQMD calculations with experimental data at RHIC energies we 
refer the reader to Refs. [321 EHl ■ 



2.3 Particle ratios — the 'horn' 

In Fig. 121 we present the excitation function of the particle ratios /n^, /n^ and (A + S°)/7r from 
central Au-|-Au (Pb-|-Pb) collisions in comparison to experimental data. The deviations between the 
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Figure 3: The excitation 
function of / , K~ / tt~ 
and (A + TP)/t[ ratios from 
5% central (AGS energies, 
SPS at 160 A-GeV and at 
RHIC energies), 7% central 
(20, 30, 40 and 80 A-GeV), 
10% central for A+E° at 160 
A-GeV Au+Au (AGS and 
RHIC) or Pb+Pb (SPS) 
collisions in comparison to 
the experimental data from 
Refs. ini 123 (AGS), 121 
123 I2S] (SPS) and ^ |2H1 
123 (RHIC) for midrapidity 
(left column) and rapidity 
integrated yields (right col- 
umn). The solid lines with 
open squares show the re- 
sults from HSD whereas the 
dashed lines with open tri- 
angles indicate the UrQMD 
calculations. The figure is 
taken from Ref. [30] . 
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transport models and the data are most pronounced for the midrapidity ratios (left column) since the 
ratios are very sensitive to actual rapidity spectra. The j-n^ ratio in UrQMD shows a maximum 
at ~ 8 A- GeV and then drops to a constant ratio of 0.11 at top SPS and RHIC energies. In the case 
of HSD a continuously rising ratio with bombarding energy is found for the midrapidity ratios which 
partly is due to a dip in the pion pseudo-rapidity distribution at RHIC energies (cf. Fig. 1 in Ref. 
|31j). The 47r ratio in HSD is roughly constant from top SPS to RHIC energies, however, larger than 
the ratio from UrQMD due to the lower amount of pion production^ and a slightly higher yield 
(cf. Fig. 12)). Nevertheless, the experimental maximum in the j-n^ ratio is missed, which we address 
dominantly to the excess of pions in the transport codes rather than to missing strangeness production. 
Qualitatively, the same arguments - due to strangeness conservation - also hold for the (A + TP)/tt 
ratio, where the pronounced experimental maxima are underestimated due to the excess of pions in the 
transport models at top AGS energies (for HSD) and above ~ 5 A-GeV (for UrQMD). Since the K~ 



^Thc lower amount of pions in HSD is essentially due to an energy-density cut of 1 GeV/fm'^, which does not allow to 
form hadrons above this critical energy density (T7j . 
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yields are well reproduced by both approaches (cf. Fig. |2)) the deviations in the K~ jn^ ratios at SPS 
and RHIC energies in UrQMD can be traced back to the excess of pions (see discussion above). 

We stress that the maximum in the (A + S°)/7r ratio is essentially due to a change from baryon 
to meson dominated dynamics with increasing bombarding energy. Similar arguments hold for the 
experimentally observed maxima in the ratio H/tt (cf. Ref. However, the 'horn' in the jn^ 

ratio at ~30 A-GeV is not described by both transport models. 

3 Transverse mass spectra — the 'step' 

We now focus on transverse mass spectra of pions and kaons/antikaons from central Au+Au (Pb+Pb) 
collisions from 2 A-GeV to 21.3 A-TeV and compare to recent data (cf. Ref. P^). Without explicit 
representation we mention that the agreement between the transport calculations and the data for pp 
and for central C+C and Si+Si is quite satisfactory no obvious traces of 'new' physics are visible. 
The situation, however, changes for central Au+Au (or Pb+Pb) collisions. Whereas at the lowest 
energy of 4 A-GeV the agreement between the transport approaches and the data is still acceptable, 
severe deviations are visible in the spectra at SPS energies of 30 and 160 A-GeV [22] • We note 
that the tt^ spectra are reasonably described at all energies while the inverse slope T of the 
transverse mass spectra in Eq. is underestimated severely by about the same amount in both 
transport approaches (within statistics). The increase of the inverse slopes in heavy-ion collisions 
with respect to collisions, which is generated by rescatterings of produced hadrons in the transport 
models, is small because the elastic meson-baryon scattering is strongly forward peaked and therefore 
gives little additional transverse momentum at midrapidity. 

The question remains whether the underestimation of the slopes in the transverse mass spectra 
P3] might be due to conventional hadronic medium effects. In fact, the my slopes of kaons and antikaons 
at SIS energies (1.5 to 2 A-GeV) were found to differ significantly j^. As argued in ^H] the different 
slopes could be traced back to repulsive kaon-nucleon potentials, which lead to a hardening of the 
spectra, and attractive antikaon-nucleon potentials, which lead to a softening of the K~ spectra. 
However, the effect of such potentials was calculated within HSD and found to be of minor importance 
at AGS and SPS energies ^H] since the meson densities are comparable to or even larger than the baryon 
densities at AGS energies and above. Additional self energy contributions stem from interactions 
with mesons; however, s-wave kaon-pion interactions are weak due to chiral symmetry arguments and 
p-wave interactions such as tt + iiT ^ K* transitions are suppressed substantially by the approximately 
'thermal' pion spectrum P7j . 

Furthermore, we have pursued the idea of Refs. [IHl ISH] that the spectra could be hardened 
by string-string interactions, which increase the effective string tension a and thus the probability to 
produce mesons at high mrp [201 EHl- In order to estimate the largest possible effect of string-string 
interactions we have assumed that for two overlapping strings the string tension o is increased by a 
factor of two, for three overlapping strings by a factor of three etc. Here the overlap of strings is defined 
geometrically assuming a transverse string radius -R^, which according to the studies in Ref. [lU] should 
be Rs < 0.25 fm. Based on these assumptions (and i?s=0.25 fm), we find only a small increase of the 
inverse slope parameters at AGS energies, where the string densities are low. At 160 A-GeV the model 
gives a significant hardening of the spectra by about 15%, which, however, is still significantly less than 
the effect observed in the data. 

Our findings are summarized in Fig. EJ where the dependence of the inverse slope parameter T (see 
Eq- ©) on is shown and compared to the experimental data [El E] for central Au+Au (Pb+Pb) 
collisions (l.h.s.) and pp reactions (r.h.s.). The upper and lower solid lines (with open circles) on the 
l.h.s. in Fig. |3] correspond to results from HSD calculations, where the upper and lower limits are due 
to fitting the slope T itself, an uncertainty in the repulsive i^^-pion potential or the possible effect of 
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Figure 4: Comparison 
of the inverse slope pa- 
rameters T for and 
K~ mesons from central 
Au+Au (Pb+Pb) collisions 
(l.h.s.) and pp reactions 
(r.h.s.) as a function of the 
invariant energy -y/i from 
HSD (upper and lower solid 
lines) and UrQMD (open 
triangles) with data from 

Refs. ini El Ha 123 123 12H1 

for AA and HH UHl EH] for 

pp collisions. The upper 
and lower solid lines result 
from different limits of 
the HSD calculations as 
discussed in the text. The 
figure is taken from Ref. 
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string overlaps. The slope parameters from pp collisions (r.h.s. in Fig. 0} are seen to increase smoothly 
with energy both in the experiment (full squares) and in the HSD calculations (full lines with open 
circles). The UrQMD results for pp collisions are shown as open triangles connected by the solid line 
and systematically lower than the slopes from HSD at all energies. 

We mention that the RQMD model jHE] gives higher inverse slope parameters for kaons at AGS and 
SPS energies than HSD and UrQMD, which essentially might be traced back to the implementation of 
effective resonances with masses above 2 GeV as well as 'color ropes' that decay isotropically in their 
rest frame A more detailed discussion of this issue is presented in Ref. |8n] . 



4 Thermodynamics in the T — /i^ plane 

This still leaves us with the question of the origin of the rapid increase of the slopes with invariant 
energy for central Au+Au collisions at AGS energies and the constant slope at SPS energies (the 'step'), 
which is missed in both transport approaches. We recall that higher transverse particle momenta either 
arise from repulsive self energies - in mean-field dynamics - or from collisions, which reduce longitudinal 
momenta in favor of transverse momenta ■ As shown above in Fig. |3] conventional hadron self- 
energy effects and hadronic binary collisions are insufficient to describe the dramatic increase of the 
slopes as a function of ^/s. This indicates additional mechanisms for the generation of the pressure 
that is observed experimentally. 

Here we propose that additional pre-hadronic/partonic degrees of freedom might be responsible for 
this effect already at ~ 5 A-GeV. Our arguments are based on a comparison of the thermodynamic 
parameters T and extracted from the transport models in the central overlap regime of Au-|-Au 
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Figure 5: The solid line charac- 
terizes the universal chemical freeze- 
out line from Cleymans et al. [TTj 
whereas the full dots (with er- 
rorbars) denote the 'experimental' 
chemical freeze-out parameters from 
Ref. jl7j. The various symbols 
stand for temperatures T and chem- 
ical potentials fiB extracted from 
UrQMD transport calculations in 
central Au+Au (Pb+Pb) collisions 
at 21.3 A-TeV, 160, 40 and 11 
A-GeV jlH] (see text). The stars in- 
dicate the tri-critical endpoints from 
lattice QCD calculations by Karsch 
et al. |1H] (left point) and Fodor and 
Katz fTI (right point). The horizon- 
tal line with errorbars is the phase 
boundary from [7]. 



collisions ^H] with the experimental systematics on chemical freeze-out configurations jl7j in the T, /i^ 
plane. The solid line in Fig. characterizes the universal chemical freeze-out line from Cleymans et 
al. jlT] whereas the full dots with errorbars denote the 'experimental' chemical freeze-out parameters 
- determined from the fits to the experimental yields - taken from Ref. jlTj. The various symbols 
(in vertical sequence) stand for temperatures T and chemical potentials hb extracted from UrQMD 
transport calculations in central Au+Au (Pb+Pb) collisions at 21.3 A-TeV, 160, 40 and 11 A-GeV jlH] 
as a function of the reaction time (from top to bottom). The open symbols denote nonequilibrium 
configurations and correspond to T parameters extracted from the transverse momentum distributions, 
whereas the full symbols denote configurations in approximate pressure equilibrium in longitudinal and 
transverse direction. 

During the nonequilibrium phase (open symbols) the transport calculations show much higher tem- 
peratures (or energy densities) than the 'experimental' chemical freeze-out configurations at all bom- 
barding energies (> 11 A- GeV). These numbers are also higher than the tri-critical endpoints extracted 
from lattice QCD calculations by Karsch et al. |1H1 and Fodor and Katz [7j (stars with horizontal 
error bars). Though the QCD lattice calculations differ substantially in the value of hb for the critical 
endpoint, the critical temperature Tc is in the range of 160 MeV in both calculations, while the energy 
density is in the order of 1 GeV/fm^ or even below. Nevertheless, this diagram shows that at RHIC 
energies one encounters more likely a cross-over between the different phases when stepping down in 
temperature during the expansion phase of the 'hot fireball'. This situation changes at lower SPS or 
AGS (as well as new GSI SIS-300) energies, where for sufficiently large chemical potentials ^b the cross 
over should change to a first order transition [4S], i.e., beyond the tri-critical point in the (T, /x^) plane. 
Nevertheless, Fig. El demonstrates that the transport calculations show temperatures (energy densities) 
well above the phase boundary (horizontal line with errorbars) in the very early phase of the collisions, 
where hadronic interactions practically yield no pressure, but pre-hadronic degrees of freedom should 
do. This argument is in line with the studies on elliptic fiow at RHIC energies, that is underestimated by 
30% at midrapidity in the HSD approach for all centralities |21j- Only strong early stage pre-hadronic 
interactions might cure this problem. 
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5 Conclusions 



Summarizing this contribution, we point out that baryon stopping jHO] and hadron production in central 
Au+Au (or Pb+Pb) colhsions is quite well described in the independent transport approaches HSD and 
UrQMD. Also the 'longitudinal' rapidity distributions of protons, pions, kaons, antikaons and hyperons 
are similar in both models and in reasonable agreement with available data. The exception are the 
pion rapidity spectra at the highest AGS energy and lower SPS energies, which are overestimated by 
both models ^7]. As a consequence the HSD and UrQMD transport approaches underestimate the 
experimental maximum of the j-n^ ratio ('horn') at ~ 20 to 30 A-GeV. However, we point out that 
the maxima in the ji:^ and (A + TP)/t: ratios dominantly reflect a change from baryon to meson 
dominated dynamics with increasing bombarding energy. 

We have found that the inverse slope parameters T for mesons from the HSD and UrQMD 
transport models are practically independent of system size from pp up to central Pb+Pb collisions and 
show only a slight increase with collision energy, but no 'step' in the transverse momentum slopes. 
The rapid increase of the inverse slope parameters of kaons for collisions of heavy nuclei (Au+Au) 
found experimentally in the AGS energy range, however, is not reproduced by both models (see Fig. HJ. 
Since the pion transverse mass spectra - which are hardly effected by collective flow - are described 
sufficiently well at all bombarding energies [201 , failure has to be attributed to a lack of pressure. We 
have argued - based on lattice QCD calculations at finite temperature and baryon chemical potential 
/i_B [3 EH] as well as the experimental systematics in the chemical freeze-out parameters |17] - that 
this additional pressure should be generated in the early phase of the collision, where the 'transverse' 
energy densities in the transport approaches are higher than the critical energy densities for a phase 
transition (or cross-over) to the QGP. The interesting finding of our analysis is, that pre-hadronic 
degrees of freedom might already play a substantial role in central Au+Au collisions at AGS energies 
above ~ 5 A-GeV . 

We recall that the systematic studies of in-plane and elliptic proton fiow in Au+Au collisions at 
AGS energies IHH 1221123 also indicate a 'softening' of the nuclear equation of state (EoS) at 4-6 A-GeV, 
which further supports our present findings. 
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